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Recently , a generalized gravity theory was proposed by Harko etal where the Lagrangian density 
is an arbitrary function of the Ricci scalar R and the trace of the stress-energy tensor T, known as 
F(R, T) gravity . In their derivation of the field equations, they have not considered conservation 
of the stress-energy tensor. In the present work , we have shown that a part of the arbitrary 
function f(R, T) can be determined if we take into account of the conservation of stress-energy 
tensor , although the form of the field equations remain similar. For homogeneous and isotropic 
model of the universe the field equations are solved and corresponding cosmological aspects has 
been discussed. Finally , we have studied the energy conditions in this modified gravity theory both 
generally and a particular case of perfect fluid with constant equation of state. 

I. INTRODUCTION 

Recent observational predictions [1] that our universe is going through a phase of accelerated 
expansion put new avenues in modern cosmology. A class of people are making attempts to accomodate 
this observational fact by choosing some exotic matter (known as dark energy) in the framework of 
general relativity. There are several choices for this exotic matter namely a) the quintessence scalar field 
models [2], the phantom field [3], K-essence [4] ,tachyon field [5], quintom [6] etc.,b) the dark energy models 
including Chaplygin gas [7] and so on. On the other hand , there are attempts to modify the gravity 
theory itself to accomodate the present accelerated phase. A natural generalization is to choose a more 
general action in which the standard Einstein - Hilbert action is replaced by an arbitrary function of 
the Ricci scalar R [8] (i.e.f(R)) and is known as f(R) -gravity . This modified theory may explain this 
late time cosmic acceleration [9]. These f(R) models can satisfy local tests and unify inflation with dark 
energy [10]. Also it is possible to explain the galactic dynamics of massive test particles in this modified 
gravity theory without any dark matter [11-15]. For detailed review of f(R) - gravity one may refer to 
[8,16]. 

Recently ,a further generalization of f(R) - gravity theory has been done by Harko etal [17]. They 
choose the Lagrangian density as arbitrary function f(R,T) where as usual R is the Ricci scalar and T 
is the trace of the energy -momentum tensor. The justification of choosing T as an argument for the 
Lagrangian is from exotic imperfect fluids or quantum effects (conformal anomaly). They have argued 
that due to the coupling of the matter and geometry ,this gravity model depends on a source term,which 
is nothing but the variation of the matter stress-energy tensor. As a result , the motion of test particles 
is not along geodesic path due to the presence of an extra force perpendicular to the four velocity . The 
cosmic acceleration in this modified f(R,T) theory results not only from geometrical contribution but 
also from the matter content. Subsequently , Houndjo [18] has chosen f(R,T) as fi(R) + fi(T) and 
discussed transition of matter dominated era to an accelerated phase. Very recently, Sharif etal [19] have 
studied thermodynamics in this f(R,T) theory and Azizi [20] have examined the possibility of wormholc 
geometry in f(R,T) gravity. 

In the present article ,we have formulated the f{R, T) gravity theory in an unorthodox manner. Though 
the action is a coupling of geometry and matter ,but still we restrict ourselves to the special cases where 
test particles move in a geodesies. As a result ,the Lagrangian has some restricted form, keeping the 
field equations same. The alternative derivation of f(R,T) gravity and some specific choice for f(R,T) 
has been presented in section II. Also admissibility of some known matter fields has been examined in 
this section . Section III deals with cosmological solutions for homogeneous and isotropic model of the 
universe with some physical interpretations. Energy conditions in this modified gravity theory has been 
examined both in a general way as well as for perfect fluid in section IV. Finally, at the end there is a 
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brief summary of the entire work in section V. 

II. f(R, T) GRAVITY THEORY : A MODIFICATION 

In this gravity theory [17], the gravitational Lagrangian density is given by an arbitrary function 
f(R, T) of two variables : One is the Ricci scalar R and the other is the trace of the energy -momentum 
tensor T(= T^g^). So the complete action of this theory is written as [17] 

A= ikiJ f{R ' T ^^~^ x + / L ^V^gd 4 x (i) 

where the stress-energy tensor of the matter (T^ u ) can be obtained from the matter Lagrangian density 
L m as [21] 



rp 2 8{^g~L m ) 

— 5g^ 



This can be simplified further assuming L. m depends only on g^ v but not on its derivatives as 

J dg^ v 



Tfiv 9[ivL m 2 (3) 



Using the standard text book result namely 

SR = R^SgT + g^DSg^ - V M V„<J<r (4) 



and the shortcut notations: 



df(R,T) d.f(R,T) 
}R - dR " /T ~ dT (5) 



the variation of the above action can be written as 



(6) 

Now performing by parts integration to the second and third terms in the r.h.s. of equation (6), one 
obtains the field equations in f(R,T) gravity theory as [17], 

f R R„ v - \f{R, T) 9llv + (g^D - V„V„)/r - 8IL7V - f T (T^ + 6^) (7) 

with 

e - =^J^ = - 2T - + ^L m - 2^gj^ (8) 

It is to be noted that if f(R,T) — f(R) then we get back to the field equations for f(R) gravity. 

Now we can proceed further ,with the field equations (7) in the following three cases: 
a)f(R,T) = R + h(T) 
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b) f(R,T) = R.h{T) 

c) f(R,T) is arbitrary 

• Case-(a) : f(R, T) = R + h(T) 
For this choice of f(R,T) the field equations (7) now simplify to 

G„„ - 8nr M , y - hr(T)(T^ + e„„) + ^/i(T) ff ^ (9) 

Now taking divergence of both sides of the above field equations (9) and assuming conservation of energy 
-momentum tensor (i.e. VT^,, = 0) we obtain 

(T^ + e^)Whr(T) + hf(T)WQ^ + |g^V^(T) = (10) 

This shows that the form of h(T) is not arbitraryit depends on the choice of the matter field. We 
consider now some known matter fields as examples: 
Example-I : Electromagnetic Field . 

The matter Lagrangian has the form 

L m = -^F a pF ia g^gP° (11) 
with F^ v ,the electromagnetic field tensor. So from equation (8) we have 

= — (12) 

As a result equation (10) simplifies to d Q^? — i.e. h(T) turns out to be a constant. Thus for this 
choice of f(R, T) electromagnetic field is not possible. 

Examplc-II : Perfect fluid. 

In case of perfect fluid, the stress-energy tensor has the usual form 

T^ v = (p + p)u^u v - pg^ (13) 

and the matter Lagrangian can be taken as L m — —p. Here p and p are the usual energy density and 
thermodynamic pressure and the four velocity satisfies i) u^u^ = 1 and ii) u^V^u^ = 0. In this case 
MJ , has the explicit form 

e^ = -2T^-pg^. (14) 
Now substituting equation (14) for Mi , into equation (10) we obtain 

(T^ + pg^)Whf(T) + hf(T)g^Wp + ^g^Wh(T) = 0. (15) 

Further, if the perfect fluid has barotropic equation of state i.e. p = ujp, co, a constant then for 
homogeneous and isotropic flat FRW model h(T) has an explicit form as (w ^ — l,±g) 

h(T) = h T a (16) 
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where a = gnj^ , h Q is an integration constant. 
• Case (b) : 

For the choice f(R,T) — Rh(T), the field equations become 

_snr ia> Rh/(T) 

h(T) h(T) ^^ + u ^>- 

After taking divergence of both sides of equation (17) and considering energy conservation relation, one 
obtains the differential equation for h(T) as 

WT ^"<m - v "'T§r >< T - + e ~> - h -wv v " e ^ - ,18) 

where Ricci scalar R is related to T and h by the relation 

R = 8UT/[hf(T)(6 + T)-h(T)}. (19) 

Note that here also for electromagnetic field h(T) is restricted by the relation 9 q^? — i.e. h(T) is a 
constant. 
• Case (c): 

Here f(R,T) is totally arbitrary except the choices in the previous two cases. We start with the 
geometric identity namely [22] 

(□V, - V„n)f(R, T) = R^Wf(R, T) 

i.e. 

V^V^V, - g^U)f{R, T) = R^WfiR, T). (20) 

Now taking covariant divergence of equation (7) and using this identity we have from conservation of 
matter field 



+ e„„)V7r + /rV8^ = 0. (21) 

Note that equation (21) is not identical in form to that of equation (10) , there is one extra term in 
equation (10). Now proceeding as before we obtain the following results : 

•I. In case of electromagnetic field equation(21) is identically satisfied and hence fx is an arbitrary 
function of R and T. Thus a general form of f(R, T) can be written as 

f(R,T) = A (R,T)+A 1 (R) (22) 
where A and A\ are arbitrary functions of arguments. 
•II. Similarly ,for perfect fluid the form of f(R,T) turns out to be 

f(R,T) = A(R)+B(T) (23) 
where A{R) is an arbitrary function of R (except A{R) = R) and B(T) has the form 

B(T) =B f exp[- f —}dT (24) 
J J P + P 
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with i? ,an integration constant. If the fluid is in barotropic nature with constant equation of state 
then we have 

B(T) = B T a , a = -J— (u; ^ -1). (25) 

Thus the choice of the function f(R,T) depends to a great extend on the matter field taken into 
account. 



III. COSMOLOGICAL SOLUTIONS AND CONSEQUENCES 

We now try to find cosmological solutions for the first choice of f(R,T) for perfect fluid in the back- 
ground of flat FRW model. The Einstein field equations are 

3H 2 = p + h Q (l-3uj) a - 1 p a (26) 

and 

2H + 3H 2 = -p+^h (l-3u) a p a (27) 
where we have used the solution (16) for h(T) . 

In Einstein gravity , the above field equations correspond to a non-interacting two-fluid system of which 
one is the usual perfect fluid (Fluid -1) that we have considered in f(R 1 T) -gravity theory while the other 
fluid system (Fluid-2) is also a perfect fluid having energy density and pressure 

Pd = h (l - 3Lo) a - 1 p a , Pd = ~ho(l - 3w) a p a (28) 



The equation of state of the additional fluid (i.e. Fluid-2) is 



with 



Pd 1 - 3w 

= — = 7T- (29) 

Pd 2 



p d +p d = y(l + 3 W )(l-3 W ) Q -y 



and 



Pd + 3 Pd = ^(9io-l)(l-3o J ) a - 1 p a , (30) 

The nature of the two non-interacting fluids in different stages of the evolution of the universe are 
shown in table I : 

Thus, although both fluids start simultaneously at the ultra- relativistic equation of state (stiff fluid) 
but fluid -2 advances more rapidly so that it reaches the quintessence equation of state when the actual 
fluid (i.e. Fluid- 1) has still positive pressure and finally fluid-2 reaches the phantom era when the fluid-1 
is in quintessence era. Further, it is to be noted that although both the fluid components have constant 
equation of state but the effective one fluid system has always variable equation of state. 
As for the normal fluid (i.e. Fluid-1) we have p = up so from the conservation of energy - momentum 
tensor i.e. 



p + 3H(p + p) = 



G 



TABLE I: Evolution of the Universe and the nature of the 2-fiuids 



Evolution of the universe 


Fluid-1 


Fluid-2 


Ultra relativistic era 
Early Universe before radiation 
Radiation era 
After radiation era 
Till dust era 
Before quintessence era 
Quintessence era 
Phantom era 


u= 1 (Stiff fluid) 
Fluid in pre-radiation era(i < u < 1) 
w= | 

5 < w < 5 
<o; < | 

-| < w < 

-1< u < -§ 

w < -1 


w d = 1 (Stiff fluid) 
Normal fluid with +ve pressure(0 < Lo d < 1) 
cjd = :dust 

- ve pressure but non-exotic in nature i.e. satisfies SEC 
Exotic fluid with — | < cud < — § 
Exotic fluid upto phantom barrier 
Phantom fluid 
Ultra phantom fluid 



we have on integration, 



P = Poa 



-3(l+«) 



(31) 



where p is an integration constant. Now substituting this value of p into the Friedmann equation (26) 
we obtain an integral equation for the scale factor 'a' as 



± (t - to) 



1 + 3^ , 

a 2 da 



/r , , 3(l-ui) , 

y/[di + d 2 a 2 



(32) 



with to an integration constant and d\ = 8I ^ Po , 
solution for 'a' with to = 0, ±1 as 



(a-1) 



. In the following we have explicit 



9d 2 
16 



(t-io) 2 



d 2 ' 



uj = 



(33) 



«o(t - to), w = 1 



and 
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di' 



w = -1. 



(34) 



(35) 



Note that on the phantom barrier ( i.e.w = —1) we have a big rip singularity at finite time t = t + 4 j^ - 
The other two solutions are the usual expanding solutions starting from the big-bang singularity at finite 
past. 

IV. ENERGY CONDITIONS IN f(R, T) GRAVITY 

The fundamental features to the singularity theorems as well as to those related to classical black hole 
thermodynamics [23] are nothing but the energy conditions which are consequences of the Raychaudhuri 
equation for expansion, namely , 



dO 1 2 



(36) 



Here 9 , and uj^ are respectively the expansion, shear , and rotation associated to the congruence 
defined by the null vector field k 1 ' and R )u/ is the usual Ricci tensor. Though the Raychaudhuri 
equation is not related to any gravity theory (it is purely a geometric statement) but it has some special 
reference to Einstein gravity. As the attractive character of gravity is reflected through the positivity 



Figure 1: Shows the graphical representation of the deceleration parameter q for the variation of ui and Q 



condition i.e. R^k^k v > (which implies that the geodesic congruences focus within a finite value 
of the parameter labeling points on the geodesies [24] ), so in Einstein gravity the above condition 
becomes T^ v h^k v > 0,which is the null energy condition (NEC). The weak energy condition (WEC) i.e. 
Tu, u v^v u > (u M ,a time-like vector) assumes the positivity of the local energy density and by continuity 
, WEC =>■ NEC . Similarly, we have two other energy conditions namely the strong energy condition 
(SEC) : (Rfj, u — \Rg i iv)v il v v > which by continuity implies NEC but not the WEC in general and 
the dominant energy condition (DEC): T flu v fJ, v L/ > and T^v" is not space-like imply locally measured 
energy density to be always positive and the energy flux is time-like or null . Also DEC =>■ WEC(&nd 
hence the NEC) but not necessarily the SEC ( For details of energy conditions see [25]) .For perfect fluid 
the above energy conditions have the explicit form: 



But difficulty arises in other gravity theories , particularly where may not be evaluated using 
the corresponding field equations. However , in a gravity theory if the Lagrangian density still have an 
Einstein - Hilbert term then it is possible to determine R^k^k". 



NEC : p+p>0 
WEC: p>0, p + p>0 
SEC: {p + 3p)>0, p+p>0 
DEC : p > Oand p±p>0 



(37) 



In the present f(R,T) gravity theory ,for the first two choices of f(R,T) (i.e. f(R,T) 
R + h(T) or Rh(T)) the field equations (9) or (17) can be written as, 




(38) 



and the energy conditions read as, 



NEC: p eff +p eff >0 

WEC: p eff >0, Peff+Peff>0 

SEC : Peff + Zpeff > 



> 



(39) 



and 



Peff + Peff > 

DEC : peff > Oand p e ff ±p e ff > , 



8 




Figure 2: Represents the variation of Zn = (1 — u)fl m + (1 + 3oj) against u and Q m . The range Zn > 

represents validity of NEC. 




Figure 3: Plots the variation of Zs = 3(1 + u>)f2 m — (1 — 3oj) against to and f2 m . The range Zs > stands for 

the validity of SEC. 



However , explicitly if we consider the perfect fluid case of the previous section then from the 
field equations (26) and (27) we have 



Pef f =p + / l0 (l-3w)("- 1 V' 

and 



Pe ff= P -\h (l-3oj) a p a (40) 



with a — 2(1+"] ■ Then the above energy conditions can be written as 
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NEC : (1 - w)O m + (1 + 3cj) > 
WEC : SameasNECan<m m > 
SEC : 3(1 + w)O m - (1 - 3w) > 
DEC : SameasWEC 

with TO = p/3H 2 , the density parameter for the matter considered in the Einstein gravity .Now the 
deceleration parameter q(= —(1 + j^j)), is related to u and Q m by the relation 

, = ( ^ )+ 3fU (1 _ w) (42) 

Thus NEC is satisfied until fluid- 1 is in the quintessence era and O m is restricted by the given inequality. 
Note that the above inequality holds for all fl m as long as the normal fluid ( i.e. fluid-l)is not exotic (i.e 
satisfies the strong energy condition) . The same is true for WEC as well as DEC. However, to satisfy 
the SEC Q m has a lower bound given by O m > . The variation of q has been plotted in figure 1 

and the inequalities for NEC and SEC are presented in figures 2 and 3 respectively. 

V. SUMMARY 

The paper deals with recently introduced f(R,T) gravity theory with the restriction of conservation 
of matter. As a result, although the form of the field equations remain same but now the test particles 
move in a geodesies and the choice of the Lagrangian function is not totally arbitrary. We have 
analyzed three possible choices for f(R,T) and examined whether two familiar matter fields namely 
electromagnetic field and perfect fluid are permissible or not in this modified gravity theory . It is found 
that electromagnetic field is not allowed in all the cases. For homogeneous and isotropic model of the 
universe , the explicit field equations are written for the modified gravity theory with f(R, T) = R + h(T) 
and it is found that the field equations are equivalent to Einstein gravity with a non-interacting 2-fluid 
system of which one is the usual perfect fluid in the modified theory while the second fluid (i.e. Fluid-2) 
is also a barotropic fluid with constant equation of state and will become exotic when the usual fluid (i.e. 
Fluid-1) is still a normal fluid. For some specific choice of the equation of state parameter of the usual 
fluid there are possible cosmological solutions of which one corresponds to big rip singularity. The graph 
(see figure 1) of q for the variation of Q m and ui shows that there is a natural transition from deceleration 
to acceleration although we have considered normal fluid (non-exotic) i.e. fluid-1 as the matter source in 
this modified gravity theory. In particular, if we consider only the baryonic matter (with VL m = .04) as 
the source of matter field then transition from deceleration to acceleration occurs when u> < .099. Thus 
with the normal fluid model in f(R,T) gravity theory, there is a natural transition from deceleration to 
acceleration as predicted by recent observations. 

Finally we have analyzed the energy conditions for the modified gravity theory in a general way. Also 
for the perfect fluid model of section - III we have shown the validity of the energy conditions both 
analytically as well as graphical. For future work, we shall deal with other aspects of the modified 
f(R,T) gravity theory. 
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